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Abstract. In several scale free graph models the asymptotic de- 
gree distribution and the characteristic exponent change when only 
a smaller set of vertices is considered. Looking at the common 
lO , properties of these models, we present sufficient conditions for the 

CN ■ almost sure existence of asymptotic degree distribution constrained 

on the set of selected vertices, and identify the chararteristic expo- 
^ ■ nent belonging to it. 
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1. Introduction 

Since the end of the nineties several complex real world networks and 
their random graph models have been investigated |11 El E] ■ Many of 
them possess the scale free property: the tail of the degree distribution 
decreases polynomially fast, that is, if Cd denotes the proportion of 
vertices of degree d, then q ~ C ■ c?"'^ holds for large values of c? p]. 7 
is called the characteristic exponent. 

If the whole network is completely known, the empirical estimator 
of the characteristic exponent may have nice properties. However, real 
world networks usually are too large and complex, hence our knowledge 
of the graph is partial. For several models of evolving random graphs 
the degree distribution and the characteristic exponent change when 
attention is restricted to a set of selected vertices that are close to the 
initial configuration [101 IHl US] ■ 

Starting from these phenomena, in this paper the degree distribution 
constrained on a set of selected vertices will be investigated, assuming 
that the graph model possesses the scale free property with characteris- 
tic exponent 7 > 1, and the number of selected vertices grows regularly 
with exponent < a < 1 . Sufficient conditions for the almost sure ex- 
istence of the local asymptotic degree distribution will be given. It will 
be shown that under these conditions the characteristic exponent of 
the constrained degree distribution is a (7 — 1) + 1. 

The proofs are based on the methods of martingale theory. Appli- 
cations of the general results to different graph models (e.g. to the 
Albert-Barabasi random tree pj) will be shown. 

In Section 2 we present the family of random graph models to be ex- 
amined and formulate the sufficient conditions. In Sections 3 and 4 we 
mention some results about martingales and slowly varying sequences 
to be applied in the proofs. Section 5 contains the proof of the main 
results, and in Section 6 we give some examples and applications. 

2. Main results 

In this section we present sufficient conditions for the almost sure 
existence of asymptotic degree distribution constrained on the set of 
selected vertices, and we describe that distribution. 

Let {Gn = (V^, -E'n))„gi^ be a sequence of evolving simple random 
graphs. Some vertices are distinghuised; let S'„ C V^ denote the set of 
selected vertices. 

We start from a finite, simple graph Gq = (Vq, Eq), this is the initial 
configuration with Vq = {ui, U2, ■ ■ ■ , ui}. Sq C Vq is arbitrarily chosen. 
For n > 1, at the nth step 

• one new vertex, t>„, is added to the graph: Vn = VoU{f 1, . . . , Vn}', 

• the new vertex gets some random edges, thus En-i C En, and 
every edge from E^ \ En-i is connected to f„; 
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• the new vertex can be added to the set of selected vertices, 
Vn G S'n is a random choice. 

The cr-field of events generated by the first n steps is denoted by J-'„. 
For V & Vn let the degree of v in Gn be denoted by deg„ {y). Fur- 
thermore, for n > 1, and (i > define 

X [n, d] = \{veVn: deg„ {v) = d}\- 

Y [n, d] = \{v eVn'. deg„ {v) = d, {v, Vn+i) E K+i}| ; 

I[n,d] = l^ if deg J.„) = c^, 
I otherwise. 

In some models it is possible that the new vertex does not get any 
edges at some steps. In other models the degree of the new vertex is 
fixed, for example, the degree of the new vertex is always 1 in random 
tree models. If the new vertex gets at least m edges at each step for 
some m > 0, then X [n, d] is at most | Vo| for all n and d < m. Thus 
we denote the minimal initial degree of the new vertex by m, and we 
consider X [n, d\ only for d > m. Of course, m = is also possible. 

2.1. Conditions on the graph model. We say that a discrete prob- 
ability distribution (a„) is exponentially decreasing if a„ < C ■ g" holds 
for all n > 1 for some C > and < g < 1 . A sequence (a„) is slowly 
varying if a[sn]/cin — )■ 1 as n — )■ oo for all s > 0. 

Throughout this paper, for two sequences (a„) , (6„) of nonnegative 
numbers, a„ ~ 6„ means that 6„ > except finitely many terms, and 
ttn/bn — > 1 as n — )■ oo. 

Now we can formulate the conditions on the graph model. 

Condition 1. X [n, rf] r^ Cd ■ n holds as n ^ oo for every d > m with 
probability 1, where (q) is a probability distribution and q is positive 
for all d > m. 

This means that asymptotic degree distribution exists in this graph 
model. Note that X [n, c?] — )■ oo as n — >■ oo almost surely. 

Condition 2. Cd ^ K ■ d~^ holds as d ^ oo for some positive numbers 
K and 7. 

This is the so called scale free property with characteristic exponent 
7. That is, the asymptotic degree distribution decays polynomially 
with exponent 7. This implies that Cd is positive for every d large 
enough, but we will need it for all d > m, this is included in Condition 

m 

Condition 3. For every n > 0, if wi,W2 G Vn and deg„ (wi) = 
deg„ (^2); then 

F{{wi,Vn+l) e En+lll'n) = P((w2,^^n+l) G En+l\ J^n) ■ 
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In other words, at each step, conditionally on the past, old vertices 
of the same degree get connected to the new vertex with the same 
probability. 

Condition 4. Y17=i I [hd] = Pd ■ n + o (n) holds as n —)> oo for every 
d > m with probability 1, where (pd) is an exponentially decreasing 
probability distribution. 

Loosely speaking, the degree of the new vertex has an exponentially 
decreasing asymptotic distribution. This trivially holds if the degree 
of the new vertex is fixed. 

Condition 5. For every d > m there exists a random variable Zd> 
with exponentially decreasing distribution such that 

P {Y [i, d]>l\J^i)<F{Zd>l), i>l,l>l. 

In many particular cases the following stronger condition is also met. 
There exists a random variable Z > with exponentially decreasing 
distribution such that 

(1) P(deg„(t;„) >/ I J-„_i) <P(Z>/), n>l,l>m. 

This is a sort of upper bound for the initial degree of the new vertex. 

Condition 6. For every d > m we have 

d 

kd = Y^ {Pj - Cj) > 0. 

j=m 

We will see later that the nonnegativity of kd follows from the previ- 
ous conditions; however, the positivity of kd cannot be omitted, as an 
example will show. 

2.2. Conditions on the set of selected vertices. Recall that Sn ^ 
Vn is the set of selected vertices in Gn- We emphasize that deg„ (f ) 
always denotes the degree of vertex v in Gn, not in Sn- 

We will need the following notations. The a-field generated by the 
first n steps and adding the edges of f„+i at the (n + l)st step is de- 
noted by J-j^. Furthermore, for n > 1 and d > m let 

X*[n,d\ = \{veSn: deg^{v) = d}\- 

Y* [n, d] = \{v e Sn ■■ deg„ (v) = d, {v, Vn+i) e ^„+i}| ; 

J 1 if w„ G Sn and deg„ (w„) = d, 
I [n, d\ = < 

I otherwise; 

T* ( ] — \^ T* \ ^1 _ J -'- if f 71 G On, 

^-^ ' 1 otherwise. 

d=m \ 

The conditions on the set of selected vertices are the following. 
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Condition 7. S'„ C S'„+i for all n > 0. 

Vertices cannot be deleted from the set of selected vertices. 

Condition 8. /* (n + 1) is J^^ -measurable for all n > 0. 

At each step we have to decide whether the new vertex is to be se- 
lected immediately after choosing its neighbours. Selecting the neigh- 
bours of a fixed vertex is an example. 

Condition 9. There exists a sequence of positive random variables (Cn) 
that are slowly varying as n -^ oo, and \Sn\ = Yll^=i ^* (0 ~ Cn ■ ■"■" for 
some a > 0, with probability 1. 

This means that the size of the set of selected vertices is regularly 
growing with exponent a > 0. 

Condition 10. For every d > m 

n n 

J2E{^[^,d]\T,^^) = {q, + o{l))J2E{I*{^)\J'.-l) 

i=l 1=1 

holds a.s. as n ^ oo, with some exponentially decreasing probability 
distribution {qd)d>m ■ 

The last condition holds if the degree of the new vertex Vn is fixed, 
or its degree and I* {n) are independent conditionally on J^n-i- In that 
case sequence qa = pd satisfies the condition. It is also possible that the 
asymptotic degree distribution of the new selected vertices is different 
from {pd) if only it decays exponentially fast. 

2.3. Description of the local degree distribution. Now we for- 
mulate the main results. 

Theorem 1. Suppose that Conditions flHTOl hold for a random graph 
model {Gn,Sn), then the limits 

X*[n,d] 
hm — — -- — = Xd 

exist for all d > m with probability 1. 

The constants Xd satisfy the following recursive equations. 

^- = c, + p^^ ' ^'^ = ^^^ {d>m + l). 

Cm Cci 

Sequence (xd) is a probability distribution, that is, it sums up to 1. 
Moreover, Xd ^ L ■ d^'' as d ^ oo with L > and 

7* = a(7-l) + l- 
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Remark 1. From the proof it is clear that with Condition\^ dropped 
the limits Xd still exist and the recursive equations remain valid. The 
role of the scale free property of the graph is just to guarantee that the 
asymptotic degree distribution constrained on the set of selected vertices 
is also polynomially decaying. 



3. Martingales 

We will extensively use the following propositions that are based on 
well-known facts of martingale theory. 

Proposition 1. Let {Mn,Qn) be a square integrable martingale with 
Ml = 0, ^0 = {0, fi}- Introduce 



A„ = ^E((M,-M,_i)'|^,_i) 



i=2 



that is, the predictable increasing process in the Doob decomposition 

of M^ . Then Mn = o ( AJ log A„ j holds almost surely on the event 

{AoQ = oo}, and Mn converges to a finite limit, as n ^ oo, almost 
surely on the event {A^o < oo} . 

This is a corollary of Propositions VII-2-3 and VII-2-4 of |14J . 

Proposition 2. Let {Mn, Qn) be a square integrable nonnegative sub- 
martingale, and 

n n 

An = EMi + Y. (^ (^^1 ^-i) - ^-i) ' ^« = 5Z^^" (^^1 ^- 



H-1, 
i=2 i=2 



,1/2 



If Bn log Bn = O {An), thcn Mn ~ An on the event {An — )■ oo}. 

This is easy to prove applying Proposition [T] to the martingale part 
of the Doob decomposition of M„. 

Proposition 3. Let Yi, 1^2, ■ ■ ■ be nonnegative, uniformly bounded ran- 
dom variables, and Qn = a {Yi, . . . ,Yn). Then the symmetric differ- 
ence of the events { ^^=i F„ < oo } and { Xl^i E {Yn\ Qn-i) < oo } 
has probability 0. Moreover, 

-^ I in ^ oo) 



Y.ZlE{Yn\Gn-l) 

holds almost everywhere on the event { Yl^=i ^n = oo } . 

This proposition follows from the Levy generalization of the Borel- 
Cantelli lemma that can be found in [13] (Corollary VII-2-6). 



6 agnes backhausz and tamas f. mori 

4. Slowly varying sequences 

In the proofs we will use the basic results of the theory of regularly 
varying sequences, see e.g. [21 El E]- 

We say that a sequence of positive numbers (/3„) is regularly varying 
with exponent fi if the following holds: 

/3n ~ 7nn^ {n -^ oo) 

where (7„) is slowly varying. 

(/3„) is regularly varying with exponent /i if and only if /3[s„]//3„ — )■ s^ 
as n ^ oo for all s > 0, see Bingham [2j. 

Proposition 4. Let {an), (/3„) be nonnegative sequences such that [an) 
is slowly varying as n ^ cxo, and n~'^/3„ — )■ 1 as n — > oo for some 
A > — 1. Then the following holds. 



y^ ftjA ~ an ^ A (n -^ oo) 



This is a consequence of the results of Bojanic and Seneta [21 [3]. 

Proposition 5. Let {an), {Pn) be nonnegative sequences such that [an) 
is regularly varying with exponent 6. 

a) Suppose ^^^^ A = Bn is regularly varying with exponent ^ > 0, 
and fi + 6 > 0. Then 

n 

V" oiiPi ~ ^- — anBn {n^ oo). 

b) Suppose Yl^=iPi — o{Bn), where {Bn) is regularly varying with 
exponent fi > 0, and fi + 6 > 0. Then 

n 

y^ aA = o {anBn) (n -^ oo) . 

Proof, a) Suppose first that 5 = 0, that is, {an) is slowly varying. 
By Bojanic and Seneta [3], for a nonnegative slowly varying sequence 
{an) there always exists another nonnegative sequence (a^) such that 
On ~ a' as n — > oo, and 



a' 



(2) lim n 1 - ^^ = 0. 

n^oo \^ a'n J 

This implies that a„+i/a„ — )■ 1 as n — ;■ oo. 
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All sequences are nonnegative, hence we have 

n n n—l i 

(3) ^ a-A = «; ^ /3j- + ^ [a[ - a.+J ^ (3^ 

i=l j=l i=l jf=l 

n-1 






as n ^ cx). 

Sequence (a^) is slowly varying. By supposition, X]j=i /^i = Tn^'^ as 
n — )■ oo with some slowly varying sequence (7n), hence i""*^ Yl^j=iPj = 
'jii^'^ as 2 — )■ oo. Since A = /x — 1 > 0, by applying Proposition H] we 
obtain that 



n f i n n ^ 

E ^ E ft ~ E ";«^- "■" ~ "Un- E ""' ~ I 

i=l j=l i=l j=l ^ 



"71 /n''' 



as n — )■ oo. Combining this with ([2]) we get that the second term on 
the right-hand side of ([3]) is o (Q;n7„n^) as n — )■ oo. 

The first term is asymptotically equal to a„7„n'^ as n — )■ oo. Thus 
we get that 

n n 

y^ ttjA ~ ^ "i A ~ an7n^^ (n -^ CX)) . 

Next, let (5 differ from 0. Let «„ = fi;„n'', and -B„ = 7„n^ with slowly 
varying sequences (k„) and (7„). We have 

n n n—l 

^,5/3. = ^^5 (B. _ B._i) = n^5„ + Y,{i'-i^ + 1)') Si 

i=l 
n-1 

i=l 

n-1 

7„n^+^ - (1 + 0(1)) (55^7/+^"' (n ^ oo) . 



«=1 i=l i=l 

n-1 

i=l 

n-1 



i=l 



(7„) is slowly varying, and A = (5 + /x — 1 > —1, thus Proposition H] 
applies, and we obtain that 



n-1 



j=l i=l 



7nn^+^ - <57n J-— = j^ Inn'-"- (n ^ oo) . 



^^^ + n 



8 AGNES BACKHAUSZ AND TAMAS F. MORI 

Let us apply the already proved particular case to (/t„) and (n'^/3„). 
Then we get that 

n n n 

1=1 i=i i=i ^ '^ 

b) We can suppose that i?„ is increasing, since for every regularly vary- 
ing sequence with positive exponent one can find another, increasing 
one, which is equivalent to it. Introduce /3^ = /3„ + -B„ — -B„-i, with 
Bo = 0. Then /?; > 0, and 

n n 

j=l 1=1 

hence it is regularly varying with exponent /i. By part a) we have 

/i 

and also 

n 

y ai {Bi - Bi^i) ~ — - — anBn- 

After subtraction we obtain that XliLi ^iPi ~ o{anBn)- D 

Proposition 6. Let ai,a2--- and 61,625- •• be nonnegative numbers 
satisfying 

1 " 
lim —y ai = K < 00, lim n6„ = 1. 

n— >-oo fl ' ^ n— J-oo 

1=1 

Then 

a) exp (X]r=i (^i^i) ^-5 regularly varying with exponent K as n ^ 00; 

b) exp (X]r=i ^l^l^i) ^-5 slowly varying as n — ?■ 00 /or every bounded 
sequence of real numbers (s„). 

Proof, a) Let rri = [tn], t > 1. We have 

m m m -. 

E «^^^= E T^+ E -[^(^ + r.)-(^-l)(i^ + r,_i)] 

■n+l i=n+l 

m m ^ m m 



t ^ — ' t 

i=n+l i=n+l i=n+l 



I ^ — ' I ^ — ' '' — ' I 

i=n+l i=n+l i=n+l i=n+l 



The first sum on the right-hand side tends to 0, since 

ai \qi\ t "^ 



(4) 



E 



i 



^ E ^~^ -~ E ai\(ii\=o{l). 

i=n+l 

The second sum is -ft' log t + (1), the third one is r^ — t„_i = o (1) 
and the last one also converges to 0. 
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b) Let m = [tn], t > 1. Now we have 



i=n+l i=n+l 

By supposition 
Qi K + iri - (i - 1) Ti^i K 



aj_ {1 + Qi) aiSi 
i 



Ti-l 



hence 



+ Tj - rj_i + 



, (1 + qi) ttiSi 
q[ = ^ (^ -^ oo) 



-)■ (z -)■ oo) , 



We can complete the proof similarly to (jl]) . 



Y.^^ 



i=n-\-l 



III, I 



i=n+l 



< 



- Yl «ik^i = o(i). □ 



i=n+l 



5. Proofs 



For sake of convenience, instead of X* [ra, rf], we consider the number 
Z* [n, d] of selected vertices with degree greater than or equal to d. That 
is, for 72 > 1 and (i > m let 

n 

(5) Z* \n,d\ = \{v G Sn : deg„ (v) > d}\ = J^X* [n,j] . 

j=d 

We also need the following notations. 

n n 

J* [n,d] = 5^r [n,j] , J[n,d] = ^/Kj] . 

j=d j=d 

First we show that Theorem [T] is implied by the following proposition. 

For all d > m we have Z* [n,d] ~ ^d |>S'n| a-s. as n — )• oo with some 
positive constants Zd- In addition. 



(6) 



Zm = l, Zd = 7 [d>m + l). 

'^d—l 



a + 



Cd-l 



It is clear that 

X* [n, d] = Z* [n, d] - Z* [n, d + 1] {n>l,d>m), 
hence 

X* [n, d] = {Zd - Zd+l) \Sn\+o{\Sn\) 

a.s. as n — 7- oo. Thus the limits 

X*[n,d] 
hm — — - — = Xd 

exist for all d > m almost surely, and Xd = Zd — Zd+i for all d > m. 

It is easy to derive the recursive equations for Xd = Zd — Zd+i from 
Zm = i and equation ([6]). The denominators are positive, because 
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Conditions dl [HI and IH] guarantee that q is nonnegative, a is positive, 
and kd is positive. 

It is also easy to check that sequence {xdj is a probabihty distribu- 
tion. We have 



and 



Xda + Xd — = Xd^i h aqd {d > m + I) . 

Cd Cd-l 



Summing up the equations above we get that 

oo oo 

^Xd = ^qd = l, 



d=m d=m 

since, by Conditions [9] and [TOl a > and the sequence (g^) is a proba- 
bihty distribution. 

The next step is solving the recursion for (x^). Set 

td = —, (^d= [[—^ — [d>m). 

It is easy to check that the recursive equations of Theorem [T] are satis- 
fied by the sequence 



d-i 

Xd 

i=0 j=i 



^ d d-1 



^d + " T^ ~~ ^j + (^ 



By Condition [21 Cd ^ K ■ d~"' holds as (i — )■ oo, and by Condition [H the 
sequence (pj) is exponentially decreasing. Hence it follows, as rf — > oo, 
that 

kd = -Yl (^j- - Pi) ~ -^ ■ zr^5 

j=m 

kd -K ■ 3^ d 

td = — " ■ ~ 



Cd K-d-' 7-1 

d-l / X d-1 






i=0 ^ "'-' i=0 



for some K' > 0. By Condition [TO] the sequence (g^) is exponentially 
decreasing, thus the series in the expression 



d 

ad (td + a) 



^^ = - /. I ^^ Yl "^^^" 
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converges. Using the asymptotics of (a^) and (t^) we get that 

1 "^ 

^d = — 7—- — r y^ aiQitt ~ L ■ (i-°(-T-i)-i 
^d [td + a) -^ 

for some L > 0. 

Consequently, the degree distribution constrained on the set of se- 
lected vertices decays polynomially, and the new characteristic expo- 
nent is determined by a and 7, namely, 7* = a (7 — 1) -|- 1, as stated. 

Therefore Theorem [1] is indeed a consequence of (j6]). 

Let us continue with the proof of ([6]). We proceed by induction on 
d. 

The case d = m is obvious, because the initial degree is never less 
than m, and the degree of a vertex cannot decrease, thus every vertex 
in Sn \ So has at least m edges. 

Suppose that 

Z*[n,d-l]r^ Zd-i\Sn\ (n -)■ 00) 

holds for some Zd-i > and d > m + 1 almost surely. 

First we determine the expected number of vertices of degree > d in 
Sn+i, given J>^, for n > 1. Every vertex in S'„ counts if its degree is at 
least d in Gn, or if its degree is equal to rf — 1 in G„ and it gets a new 
edge from f„+i. The new vertex Vn+i counts if it falls into Sn+i and 
its degree is > d in Gn+i- Thus the following equality holds for every 
n>l. 

(7) Z* [n + l,d] = Z* [n, d] + Y* [n, d - I] + J* [n + l,d] . 
Taking conditional expectations with respect to J-'„ we obtain that 

(8) E{Z*[n + l,d]\J^^) 

= Z* [n, d] + E{Y*[n,d- 1]| J-„) + E{J*[n + l, d\\J^n) • 

By Condition [3l vertices of the same degree are connected to f „+i with 
the same conditional probability. This implies that 

X [n, d] may be equal to zero, then Y [n, d] = as well. We will consider 
all quotients of the form 0/0 as 1. 

The middle term on the right-hand side of (jHj) can be transformed 
by the help of (^. 

(10) E{Z*[n + l,d]\J^n) 

= Z* [n,d]+X* [n,ci- 1] ^^^J^'^r \J'-^"^ +E( J* [n+l,d\\J-r.) ■ 

X [n,d — 1\ 
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By (E]), X* [n, d- 1] = Z*[n,d - 1] - Z* [n, d\, hence from equation 
( TTU]1 we obtain that 

(11) E{Z* [n + l,d]\:F^) = Z* [n,d] ^1 - ^ ^^ [^ V-^}j'^"0 

+ z-Krf-i] ^^yj|;^^~_^||'^"^ +E(j*[n + i,rf]|j-„) 

for all n > 1. 
For z > 1 define 



if X [i, d] = 0; 

X[i,c[\ 



i-^'^IF) ' if.f|M]>o. 



Set c [1 , (i] = 1 and for n > 2 define 



n-l 



12) c[n,d] = Y[b[i,d]. 



i=l 



Then for n large enough we have 



c[n,d] [^ E{Y[n,d]\J^^ 



c[n + l,(i] \ X[n,d] 

For several particular models it is quite easy to compute the con- 
ditional expectations E {Y [i,d — 1]| J-^), and hence, to determine the 
asymptotics of c[n, rf]. In the present general case the conditional ex- 
pectation is not specified. However, as the following sequence of lem- 
mas shows, the asymptotics of the partial sums can be described, and 
one can calculate the asymptotics of c [n, d\. The proof of the lemmas 
will be postponed to the second part of this section. We emphasize 
that in the lemmas the induction hypothesis is assumed all along. 

Consider the partial sums 

n 

S[n,d] = Y,E{Y[i,d]\J',) {n>l). 

i=l 

Lemma 1. For all d > m we have 

n 

(13) S[n,d] = J2E{Y[i,d]\J^i) = kd-n + o{n) {n ^ oo) 



i=l 



with probability 1. 



Remark 2. It is clear from the definition that S[n,d] is nonnegative, 
hence LemmaUl immediately implies kfi> for alld > m (cf. Condition 

Lemma 2. 

c [n, d]r^ a [n, d] ■ n^<i/''<i (n -> cx)) , 
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a.s. for all d > m, where a[n,d] is positive and slowly varying as 
n — )■ oo. 



By equation flTTl) . the process 

(14) V[n,d] = c[n,d-l]Z*[n,d] {n > 1) 

is a submartingale. Let A [n, d] denote the increasing process in the 
Doob decomposition of V [n, d] ; it is given by 

(15) A[n,d] = J2c[^ + l,d-l] Z* [t, d - 1] ^ ^y|: "I^I^jj ^'^ 

n 

+ J2c[i + l,d-l]E{r[i+l,d]\J^,). 



i=l 



i=l 



First we describe the asymptotics of A [n, d]. 

Lemma 3. Suppose that Z* [n,d — 1] ~ za-i \Sn\ holds a.s. for some 
d> m + 1, as n ^ oo, then 



kd-i 

A [n, d] ~ '^ ^"''' ' fc^7'^"^' ^ [^' ^] Cn^"^^'-'^''-' a.s. 






Next, we compute an upper bound for the conditional variances. 
Define 

n 

5[n,d] = ^Var(V[z,rf]| J-i^i) {n > 2) . 

i=2 

Lemma 4. Suppose that Z* [n,d — 1] ~ z^-i \Sn\ holds a.s. for some 

1 /9 

d > m + 1, as n ^ oo, then B [n, d] ' log 5 [n, d] = O {A [n, d\). 

Therefore Proposition [2] implies that V [n, d] ~ A [n, d] almost surely 
as n — )■ cxD. Finally, by Lemma [2] and Lemma [3] we obtain the asymp- 
totics 

Z* [n, d] ""'-' . / Cnn" (n ^ oo) . 

a + -^^ 

Cd-l 

Consequently, we have 

Z* [n, d] ~ ZdCnn"" {n -^ cx)) 
with 

kd—-\ I V^CXD 

16 Zd = r— • 

Cd-l 

The size of Sn is asymptotically equal to Cn'^" by Condition |9l Thus 
the proof of ([6]) can be completed by using Lemmas [THU D 

Now we continue with the proofs of Lemmas [THU 
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Proof of Lemma [T]. Similarly to equation ([7]), but considering all 
vertices, we see that 

X [i + 1, j] = X [z, j] - Y [z, j] + Y[i,3-l]+I[i + 1, j] 

for every i > and j > m. Adding up for z = 1, . . . , n we obtain that 

n n n+1 

(17) X[n+l,J]-X[l,J] = -Y,y[hJ] + Y.^[^,J-l] + J2^l'^3] 

for every j > m and ra > 1. By Conditions [1] and IH from f lT7|) it follows 
that 

n n 

X] ^ [^'^ ~ 1] ~ X^ ^ [^'^1 = (^i - Pi) ■ ^ + o H 

holds almost surely, as n — )■ oo, for every j > m. Adding this up for 
j = m, . . . ,d we get 

n d 

(18) 22 ^ [^' "^1 = ~ z2 ^^^ - Pj) ■ n + o (n) = kd ■ n + o (n) 

«= 1 j=m 

a.s., as n — 7- oo. Therefore it is sufficient to prove that 
1 " 

-^(y[i,rf]-E(r[i,d]|j^,)) ^0 (n^oo). 
^ i=l 

Fix d> m, and for n > 1 let M„ = J2Li (^ [^, ^] - ^ (^ [^, d]\J^^)), 
Gn = -^n+i- It is clear that (M„, ^„) is a martingale. Using Condition^] 
we will derive an upper bound for the corresponding increasing process 
An introduced in Proposition [1] 

n n 

(19) A„ = ^Var(r[2,d]|J-,)<^^(r[2,rf]'|j;,) < 

n n 

i=l i=l 

n n 

<Y,CrE{Y[i,d]\F.;) + Y,E{Zll{Za>C,)) 

i=l i=l 

for any Ci > 0. Fix e > such that k = E (e'^^'') is finite, and for i >3 
choose Ci = -log 2. The function z h- >• z^e~^^ is decreasing for z > -, 
hence z'^e^'^^ < Cf e^^'"* for z > Ci. This implies 

(20) E [Zll {Zd > C,)) < C^-^^'E [e'^'I {Z, > C,)) 

The infinite sum of these terms converges, thus the second sum on the 
right-hand side of (fT9|) is bounded for fixed d. 
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On the other hand, F [i, c?] < \Vi\ <i + l follows from the definition, 
therefore 

J2C,E{Y[z,d]\T,)<Y,-^ogf{t + l) = 0{nHogn). 

i=3 i=3 

Thus An = O (n^ logra). This bound can be further improved as follows. 
Applying Proposition [1] to the martingale (M^) we get that M„ = 
O (77,^+'') a.s. for all i] > 0. Equation flTSjl imphes that Y17=i ^ [^' ^] ~ 
O (n), therefore 



4 = 1 4 = 1 

= Cj f^F [«,f/] - mA = O (n^+^logn) . 

We obtain that An = O [n^^'^logn). Hence by Proposition [1] we have 

Mn = o (n2^''^logn] a.e. on the event {Aoo = 00}, for all i] > 0. 

Therefore M„ = o{n) holds almost surely, and this completes the proof 
of Lemma [H D 



Proof of Lemma [21 Fix an arbitrary d > m. Lemma [T] and the 
induction hypothesis imply that 

E{Y[n,d]\J^n) S[n,d]-S[n-l,d] ^ . , 

vr ^1 = v\ A^ ^ [n^ 00). 

X [n, d\ X [n, d\ 

Thus in flT2l) we can apply the approximation 1 — x = e~^~^ ^^ ^ (x — )■ 
0). Set 

^ _ S[i,d]- S[i-l,d] _ Cd 



Cd X [z, d] 

a X [i,d] > 0, and ai = bi = otherwise. Then a^ and bi are nonnega- 
tive. In addition, 

}Lya, = —y{S[z,d]-S[z-l,d]) = ^^!^^'^ 
n^ CdU^ CdU Cd 

4 = 1 4=1 

as n — )■ 00, by Lemma [TJ According to the induction hypothesis, 
X [n,d] ^ Cd ■ n, which implies that n6„ — )■ 1 as n — )■ 00. Therefore 
Proposition E] applies to the sequences (a„) and (6„) with K = kd/cd- 
Thus, due to part a). 



n \ / n 

exp 



V- ,\ ^ S[t,d]-S[i-l,d] 



.4 = 1 / \4=1 

is regularly varying with exponent K. 
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The remainder terms produce a slowly varying function, because by 
part b) of Proposition we get that 

(^ 2,2 ^ (^ f S[i,d] -S[i- l,d] y 
exp [l.-^b,^^j = -P [1. [ YM J "' 

is slowly varying supposed the sequence (sj) is bounded. 

From these the asymptotics of c [n, d] readily follows. D 

Proof of Lemma [31 By f lT^ . A[n, d] = Ai + A2, where 

n 

A2 = ^c[i + l,d - l]E {J* [i + l,d]\7^) . 

i=l 

Here we already know the asymptotics of S[n,d — 1] and c[n,d — 1] 
from Lemmas [Hand [21 In addition, Z* [n, d — 1] ~ Zd-i |5'„| ~ Zd-iCnn°' 
a.s., due to the induction hypothesis and Condition [91 It is clear from 
the definition that Y [n,d — 1] is nonnegative, thus we have 

n 

A~ Va[z + 1,^-1] z'=^-i/^^-^rf_iO^" :E{Y[t,d-l]\T,) 

n 

= ^^ V a [z + 1, ci - 1] 1^^-^/-^-^+^-^ E{Y[i,d-l]\J^,). 
1=1 

Let us apply part a) of Proposition [5] in the following setting. 

an = a[n + l,d-l]Cn, Pn = E{Y[n,d- l]\J^n) 
for n > 1; 

= \- a — 1. 

Cd-l 

Condition M and Lemma [2] guarantee that q;„ is slowly varying. Fur- 
thermore, Yl^=i Pi — S[n,d — 1] ~ kd-iu as n — )■ 00 with k^-i > 0, 
hence 7„ = kd-i and /i = 1 satisfy the conditions. Finally, fi + 6 = 
kd-i/cd-i + a > 0, because Cd-i, kd-i and a are positive due to Con- 
ditions [H [6l and[9l 

Applying Proposition [5] we obtain that 

A, ~ ^^^^i^ ■ -^^ .a[n,d-l] Cn n^-^/^-^+" 

almost surely as n — )■ 00, where z^-i, fcd_i and Cd-i are positive. 
Now we examine the second term in A[n, d\. Since 



j=m 
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we have 

d-i 

Hence by Lemma [2] 

n 

^2 ~ ^ a [z + 1, d - 1] i'^'^-i/^'^-i : 



X 

i=l 

d-l 



j=m 



Set an = a [n + 1, d - 1], 6 = kd-i/cd-i, and /3„ = E (/* {n + 1)| -Fn). 
By Proposition |3] and Condition |9] we have 

n n 

(21) J] A ~ 5^r (^ + 1) = l^n+ll - /* (1) ~ Cn^i" (n^Oo). 

i=l i=l 

Thus we can apply part a) of Proposition [5] with /i = a > 0. Assump- 
tion (5 + /i > is satisfied. Therefore we get that 

n 

y"c[t + l,d-l]r{t + l) %-- ■a\n + l,d-l]Cn n^+'^^-i/'^'^-i 

7^ a + ^^^^ 

almost surely as n — )■ oo. 

On the other hand, for a fixed j < rf — 1 we have 

n 

^c[^ + i,rf-i]i?(r[^ + i,j]|j-,) 

~ ^a [i + 1, d - 1] z^^^-i/^d-i^ (/* [i + 1, j]| j:,) 



j=l 



i=l 



by Lemma [2l In this case a„ remains the same as before, and we set 
f3n = E {I*[n + l,j]\J^n)- Using Condition [10] and equation (^T^ we 
obtain that 

n n 

n 

(g, + o (1)) J] E (r (z + 1) I J-.) = (g, + o (1)) Cnn" 



j=l i=l 



i=l 



almost surely as n — )■ oo. Thus we can apply part a) or part b) of 
Proposition [5] with /i = a, according that g^ vanishes or it is positive. 
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Then we get that 



■i 1 


i + l,d- 


-i]E{r 


[i + hj] ^i) 








_ aqj + (1) 

a + ^"'^ 



a[n,d- 1] Cnn'^+^''-^/^-^- 

Cd-l 



almost surely as n — )■ oo. Hence we conclude that 



a 



"+fcd-i/cd- 



(22) ^2- 1- Vg, + o(l) a[n,d~l]Cnr^ 

almost surely as n — )■ oo. Since (g^) is a probability distribution by 
Condition [T0| it follows that (l — X]i=m^i) ~ J^'jLdlj- This com- 
pletes the proof. D 

Proof of Lemma (4]. From equation flT^ it follows that 

n n 

E [n, rf] = ^ Var ( \/ [z, d] | J'i-i) = ^ c [z, rf - 1]^ Var ( Z* [i, rf] | Ti^i) . 

i=2 i=2 

By equation ([7]) we have 

Var ( Z* [t, d] I J-i_i) < E ( (Z* [2, d] -Z*[i- 1, d])' I J-i_i) 
= E{{Y*[i-l,d-l] + r [t, d]f\ j-i_i) 
< 2E (r* [z - 1, rf - 1]' I J-i_i) + 2E ( J* [z, d] I J-,_i) . 

Hence 

n 

(23) 5[n,t/] <2^c[z,rf-l]2E(F*[i-l,t/-l]^| J-i^i) 

=2 

n 

+ 2 ^ c[z, ci - 1]' E ( J* [2, rf] I j:,_i) = 2Bi + 252. 



i=2 



j=2 



We will estimate -Bi and B2 separately. 

Similarly to the proof of Lemma [H fix a positive e > such that 
K = E [e^^'^-^^ < 00, and set Ci = -logz. Using Condition [5] and 
inequality Y* [i,d] < Y [i,d] one can see that 

E{Y*[t-l,d-lf\j^,^,) <QE{Y*[t-l,d-l]\Ti^i) 

+ E{Zl^J{Za-i>Q)) 
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holds. For estimating the first term on the right-hand side we make 
use of equation (^ . 



E{Y*[t-l,d-l]\J',_,) 



X[i-l,d-l] 



X*[i-l,d- 1] 



E{Y[^-l^d-l]\T^ 

- x[t-i,d^i] '''*-'' 

To the second term we can apply (I2U]) : it is O ((logi)^z~^). 
From all these we obtain that 

Bi<2^c[t,d-l\ C, |S,_i| 



X[i- l,d-l] 



i=2 

/ n 

t;[f, u — ±j i^iog'f 

.1=2 



Note that the second sum is convergent here. In the first sum c[z, d—\\ 
can be estimated by Lemma [21 l^j-il by Condition [HI and X[i — 1, d — 1] 
by Condition [TJ In this way we obtain that 

is regularly varying with exponent b = 2kd-i/cd-i + a — l. On the other 
hand, by Lemma [T] the sum of E {Y[i — l,d — 1]\ J^t-i) is regularly 
varying with exponent 1. Therefore part a) of Proposition O implies 
that 

Bi = {a[n,d-lf{\ognfCnn''^^^''-^'''-'-^) . 
For B2 let us apply part a) of Proposition [5] with ctj = c[z, d — 1] 
and Pi = c[i, d — 1]E {J*[i, d]\ J-i-i). The regular variation of J2 Pi has 
already been proven in fl22|) . Thus, 

Returning to fl23l) we conclude that 

B [n, d]=0 (n"+2fcd-i/cd_i+r,^) 

for all ?7 > 0. Consequently, 

B [n, df'^ log B [n, d] = (n"/2+'^d-i/<=^-i+^) {n ^ 00) . 
Now the proof can be completed by comparing this with Lemma 

m □ 

Remark 3. Since S [n, d] is clearly nonnegative, LemmaUl implies that 
^d > for all d > m. This means that 

d d 

^Cj>^Pj {d>m). 
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Loosely speaking, the degree of a typical vertex is asymptotically larger 
than or equal to the degree of the new vertex. This is in accordance with 
the fact that the degree of a fixed vertex cannot decrease. 

Similarly to Lemma 1, one can prove that 

d d 

^Xj>Y^ qj {d>m), 

j=0 j=0 

which means the same for the selected vertices. 

6. Graph models 

In this section we briefly review some scale free random graph models 
and sets of selected vertices to which the results of the previous section 
can be applied. 

6.1. Generalized plane oriented recursive tree. We start from 
one edge, and at each step one new vertex and one new edge are added 
to the graph. At the nth step the probability that a given vertex of 
degree d is connected to f„ is (rf + /3) /T„_i, where f3 > —1 is the 
parameter of the model, and T„_i = (2 + /3) (n + 1) + /3. These kind of 
random trees are widely examined, see for example [3,[T5]. f3 = gives 
the Albert-Barabasi tree pl^. 

We fix an integer j > 1. At the nth step f„ is added to the set of 
selected vertices if it is at distance j from ui in Gn- Thus Sn is the jth 
level of the tree Gn- 

It is well known [9] that Condition [1] is satisfied with 

(2 + /3)r(c^ + /3)r(3 + 2/3) 



Consequently, 

Cd 



(2 + /3)r(3 + 2/3) _(3,,) 

r (1 + /3) Id ^ ooj 



and 7 = 3 + /3 satisfies Condition [21 It is clear from the definition that 
Condition [3] holds, and since the degree of the new vertex is always 
1, we have m = 1, and conditions HI O and [10] are trivially satisfied. 
Using that pd = ior d j^ 1 and pi = 1, Condition [6] is also easy to 
check. 

The distance of Vn and ui does not change after generating the edges 
from Vn at the nth step. This guarantees Conditions [7] and [81 The 
results of [TU] show that Condition [HI is satisfied with a = 1/ (2 + (3). 
It is proven that 

\Sn\- nC ^._ ., e-^(") X n^^ (logn)^'"^ (n ^ oo) , 
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where C is a positive random variable, and /x (n) = |i| log n [101 The- 
orem 2.1]. 

Thus Theorem [1] applies: the asymptotic degree distribution con- 
strained on a fixed level of the tree does exist. The new characteristic 
exponent is the following (cf. [TO', Theorem 3.1]). 

^* = a(^-l) + l = _l_(3 + /3-l) + l = 2. 

6.2. Independent edges. We start from one edge. At the nth step, 
independently of each other, every old vertex is connected to the new 
one with probability Xd/Tn-i, where d is the degree of the old vertex 
in Gn~i, < A < 2 is a fixed parameter, and T„_i denotes the sum of 
degrees in Gn-i- The restriction on A guarantees that the probability 
given above belongs to [0, 1]. It is clear that m = 0. 

We fix one vertex, v, and Sn consists of its neighbours in Gn- 
In [HI Theorem 3.1.] it is proven that the asymptotic degree distri- 
bution is given by 

2 '^ 



where 



A'^ 



Clearly, q ~ 2A (2 + A) d'^ {d — )■ oo). Thus the first two conditions 
are satisfied, and 7 = 3. Condition [3] holds, because the probability 
that a given vertex gets a new edge depends only on its actual degree. 
It is also clear that Conditions [7] and [H] hold. Condition IH] is a corollary 
of [m Theorem 2.1], and we have a = 1/2. 

In this case the initial degree of the new vertex is not fixed. It is 
proven in [8] that 



E 



\E{I[n + l,d]\Tn)-Pd\ ->0 

k=0 

almost surely as n — )■ oo. This, and the fact that {pdj is a Poisson 
distribution with parameter A imply Condition 21 

Note that the conditional distribution of Y [n, d] is binomial of order 
X [n, d] and parameter Xd/Tn < 1. One can check Condition El with Z^ 
having a suitable Poisson distribution. 

Condition [TUl can be verified basing on the fact that the degree dis- 
tribution of a new selected vertex is similar to the distribution of a 
new vertex because of the independent random choices, and the fol- 
lowing results. Theorem 2.1 in [12] states that T„ = 2An -|- o{n^~^) 
almost surely if e > is sufficiently small. Moreover, Theorem 2.2 
there implies that the maximum degree after n steps is O (i/n) almost 
surely. 
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Our Theorem [T] can be applied, so the almost sure asymptotic degree 
distribution constrained on the neighbours of a fixed vertex exists. The 
new characteristic exponent is given by 

7* = a(7-l) + l = ^(3-l) + l = 2 

(cf. dH Theorem 3.1]). 

Let us modify this example in such a way that vertices of degree 1 
never get new edges. Let 

n 

Tn-1 = ^X[n- l,d\d, 

d=2 

and choose Sn to contain all vertices of degree L Then we can see 
that the all conditions hold except Condition [HI but a;^ = for d > 1. 
This shows that positivity of k^ cannot be relaxed in order to obtain a 
polynomially decreasing degree distribution. 

6.3. Random multitrees. For M > 2 an M-multicherry is a hyper- 
graph on M + 1 vertices. One of them, called center, is distinguished, 
it is connected to all other vertices with ordinary edges (2-hyperedges), 
and the remaining M vertices form an M-hyperedge, called the base. 

We start from the complete graph of M vertices; the vertices form 
a base. Then at each step we add a new vertex and an M-multicherry 
with the new vertex in its center. We select the base of the new multi- 
cherry from the existing bases uniformly. Finally, we add M new bases 
by replacing a vertex in the selected base with the new center in all 
possible ways. 

The degree of the new vertex is always M , thus m = M. 

Let Sn be the set of vertices that are at distance j from the initial 
configuration. 

It is shown in [13] that Conditions [U [2|, and M are satisfied with 



7 = 2 + jg3Y and a = ^^^^- The other conditions are easy to check, 
using that distances in the multitree do not change. 
Therefore Theorem [1] applies, and 

M-l f 1 \ 

Another option for the set of selected vertices is the following. Fix 
an integer 1 < k < M and k different vertices. Let Sn be the set of 
vertices that are connected to all of them. Since the model is the same, 
we only have to check the conditions on the set of selected vertices. 
Now Conditions [TJ [H and [10] clearly hold. Condition [9] can be proven 
by slight modifications of the proofs of [13]. In this case 7 = 2 + j^zii 
a = 1 — -g, and 
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7. CONCLUSIONS 

We presented sufficient conditions for the existence of the asymptotic 
degree distribution constrained on the set of selected vertices. Scale free 
property and regular variation of the size of the set of selected vertices 
were essential. The new characteristic exponent depended only on 7 
and a. 

We reviewed several models satisfying these conditions and identi- 
fied their characteristic exponents applying our main result. In these 
models 7* < 2 and 7* < 7, thus the characteristic exponent decreased. 
One reason for that is the following. The selected vertices are closer to 
the initial configuration in some sense. There are more "old" vertices 
among them and their degree is larger than that of the "typical" ones. 
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